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Abstract. Motivated by structured parasite populations in aquaculture we consider a class of 
size-structured population models, where individuals may be recruited into the population with 
distributed states at birth. The mathematical model which describes the evolution of such a pop- 
ulation is a first-order nonlinear partial integro-differential equation of hyperbolic type. First, we 
use positive perturbation arguments and utilise results from the spectral theory of semigroups to 
establish conditions for the existence of a positive equilibrium solution of our model. Then, we 
formulate conditions that guarantee that the linearised system is governed by a positive quasicon- 
traction semigroup on the biologically relevant state space. We also show that the governing linear 
semigroup is eventually compact, hence growth properties of the semigroup are determined by the 
spectrum of its generator. In the case of a separable fertility function, we deduce a characteristic 
equation, and investigate the stability of equilibrium solutions in the general case using positive 
perturbation arguments. 
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1. Introduction 



In this paper, we study the following partial integro-differential equation 



d_ 

Of 



■p(s,t) + — {y{s,P{t))p{s,t)) 
j(0,P(t))p(0,t) 



p{s,P{t))p{s,t) + 



Jo 



/ P(s,y,P(t))p(y,t)dy, (1.1) 



= 0, 



(1.2) 



p(a,0)=po(a), P(t) = / p{s,t)ds 



(1.3) 



Here the function p = p(s, t) denotes the density of individuals of size (or other developmental 
stage) s at time t with m being the finite maximal size any individual may reach in its lifetime. 
Vital rates fj, > and 7 > denote the mortality and growth rates of individuals, respectively, and 
both depend on both size s and on the total population size P(t). It is assumed that individuals 
may have different sizes at birth and therefore ft(s, y, ■ ) denotes the rate at which individuals of 
size y give rise to individuals of size s. The non-local integral term in ( 11.11) represents reproduction 
of the population without external driving of the population through immigration. We make the 
following general assumptions on the model ingredients 



Our motivation to investigate model ( |1.1I) - (I1.3I) is the modelling of structured parasite popu- 
lations in aquaculture. In particular we are interested in parasites of farmed and wild salmonid 
fish that have particular relevance both industrially and commercially to the UK. These species 
are subject to parasitism from a number of copepod (crustacean) parasites of the family Caligi- 
dae. These sea louse parasites are well studied with a large literature: below we draw attention to 
some recent key review papers. Sea lice cause reduced growth and appetite, wounding, and sus- 
ceptibility to secondary infections 0, resulting in significant damage to crops and therefore they 
are economically important. For salmon, louse burden in excess of 0.1 lice per gram of fish can 
be considered pathogenic 0. The best studied species is Lepeophtheirus salmonis, principally a 
parasite of salmonids and frequent parasite on British Atlantic salmon (Salmo salar) farms [|22l . It 
also infects sea trout {Salmo trutta) and rainbow trout (Oncorhynchus mykiss). The life history of 
the parasite is direct, with no requirement for intermediate hosts. It involves a succession of ten 
distinct developmental stages, separated by moults, from egg to adult. Initial naupliar and cope- 
podid stages are free living and planktonic. Following attachment of the infectious copepodid to 
a host, the parasite passes through four chalimus stages that are firmly attached to the host, before 
entering sexually dimorphic pre-adult and adult stages where the parasite can once again move 
over the host surface and transfer to new hosts. 

The state of the art for population-level modelling of L. salmonis is represented by Revie et 
al. Il20"l . These authors presented a series of delay-differential equations to model different life- 
history stages and parameterised the model using data collected at Scottish salmon farms. A similar 
compartmental model was proposed by Tucker et al. [|2TI . The emphasis of these papers was not 



H E C^^m] x [0,oo)), p E (^([O.m] x [0,m] x [0,oo)) 
P, /i>0, 7 G C 1 ([0,m] x [0,oo)), 7>0. 



(1.4) 
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however, in analytical study, but on numerical simulation and parameterisation using field EOl and 
laboratory ETI data. An earlier model by Heuch & Mo |fT3l investigated the infectivity, in term of 
L. salmonis egg production, posed by the Norwegian salmon industry, using a simple deterministic 
model. Other authors have considered the potential for long-distance dispersal of mobile parasite 
stages through sea currents lfT8l . looking at Loch Shieldaig in NW Scotland, a long-term study site 
for sea louse research. 

In this paper, we focus on the dynamics of individuals at the chalimus to adult stages. Though 
individuals pass through a series of discrete growth stages by moulting, this outward punctuated 
growth disguises a physiologically more smooth growth process in terms of the accumulation of 
energy, and by 'size' in this paper we presume accumulation of energy, rather than physical dimen- 
sion. Sea lice reproduce sexually; however at the chalimus stage individuals are not yet sexually 
differentiated. Fertility rates thus must be considered as applying to the population as a whole, 
rather than as is usually the case the female fraction of the population. Individuals entering the 
first chalimus stage from the non-feeding planktonic stages are distributed over different sizes, 
hence we have the zero influx boundary condition (11.21) and the recruitment term in (11.11) . Our 
aim here is to present a preliminary step towards the analysis of the more complex problem of 
modelling the whole life cycle of sea lice by giving a mathematical treatment of a quite general 
scramble competition model with distributed states-at-birth. We use the term scramble competition 
to describe the scenario where individuals have equal chance when competing for resources such 
as food (see e.g. JH). Therefore all vital rates, i.e. growth, fertility and mortality depend on the 
total population size of competitors. In other populations, such as a tree population or a cannibalis- 
tic population, there may be a natural hierarchy among individuals of different sizes, which results 
in mathematical models incorporating infinite-dimensional nonlinearities, see e.g. |[T0l [TT1 . The 
analysis presented in this paper could be extended to these type of models and also to other models 
such as those that involve a different type of recruitment term. 

Here, we consider the asymptotic behaviour of solutions of model (ll.ll) - (ll.3l) . Our analysis 
is based on linearisation around equilibrium solutions (see e.g. lfT0l[T9l ) and utilises well-known 
results from linear operator theory that can be found for example in the excellent books QlElUll. 
We also utilise some novel ideas on positive perturbations of linear operators. For basic concepts 
and results from the theory of structured population dynamics we refer the interested reader to 

Bids 

Traditionally, structured population models have been formulated as partial differential equa- 
tions for population densities. However, the recent unified approach of Diekmann et ah, making 
use of the rich theory of delay and integral equations, has been resulted in significant advances. 
The Principle of Linearised Stability has been proven in [O [8J for a wide class of physiologically 
structured population models formulated as delay equations (or abstract integral equations). It is 
not clear yet whether the models formulated in [8l as delay equations are equivalent to those 
formulated as partial differential equations. 

In the remarkable paper [3], Calsina and Saldana studied the well-posedness of a very general 
size-structured model with distributed states-at-birth. They established the global existence and 
uniqueness of solutions utilising results from the theory of nonlinear evolution equations. Model 
(ll.ll) - (ll.3l) is a special case of the general model treated in [3], however, in qualitative questions 
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were not addressed. In contrast to (3), our paper focuses on the existence and local asymptotic 
stability of equilibrium solutions of system ( |1.1|) -( |1.3I) with particular regards to the effects of 
distributed states-at-birfh compared to more simple models we addressed previously, e.g. in IITOll . 
First, we establish conditions in Theorem [6] that guarantee the existence of equilibrium solutions, 
in general. Then, we show in Theorem [8] that a positive quasicontraction semigroup describes 
the evolution of solutions of the system linearised at an equilibrium solution. Next, we establish 
a further regularity property in Theorem [121 for the governing linear semigroup, which allows 
one to investigate the stability of positive equilibrium solutions of (|1.1|) - (|1.3I) . We use rank-one 
perturbations of the general recruitment term to arrive at stability/instability conditions for the 
equilibria. Finally we briefly discuss the positivity of the governing linear semigroup. 



2. Existence of equilibrium solutions 

Model (ll.ll) - (ll.3l) admits the trivial solution. If we look for positive time-independent solutions of 
( |1.1|) -( |1.3I) we arrive at the following integro-differential equation 



7(s,P*K(s) + (7 s (s,P*) + fi(s,P*))p*(s) = / (3{s,y,P*)p*{y)dy (2.1) 

Jo 

7 (0,P*)p*(0) = 0, P*= / p*(s)ds. (2.2) 

Jo 

2.1. Separable fertility function 

In the special case of 

P(s,y,P) = p 1 (a,P)p 2 (y), s,ye[0,m], Pe(0,oo), (2.3) 

where the distribution of offspring sizes is dependent upon the level of competition P, but the 
mature size at which individuals reproduce is not, equation (12.11) reduces to 

7 (s, P*K(s) + (%(s, P*) + P.))p*{s) = Pi{s, P*)P*, (2.4) 

where 

/■m 

P* = / /3 2 (y)p*(y)dy. 

Jo 

The solution of (12.41) satisfying the initial condition in (12.21) is readily obtained as 



where 
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Multiplying equation (12.51 ) by {3 2 an d integrating from to m yields the following necessary con- 
dition for the existence of a positive equilibrium solution 

1 = [ m p 2 (s)F(s,P*) f S P /i V ^* ] p , dyds. (2.6) 

Jo Jo F{y,P*n{y,P*) 

Therefore we define a net reproduction function R as follows 

fi( P) = r fftfe^uj- d, d , (2.7) 



JO 



1(S,P) 1 I Jy 1(Z,P) 



It is straightforward to show that for every positive value P* for which P(P*) = 1 holds, formula 
(12.51) yields a unique positive stationary solution where P* may be determined from equation 
(1231) as 

- P* 
Then it is straightforward to establish the following result. 

Proposition 1. Assume that the fertility function (3 satisfies (12.31) and that the following conditions 
hold true 

P{s,y,0)> fi{s,0), s,y E [0, m], P E (0, oo); J exp |- jf dyj ds < m - 1, 

(2.8) 

/ /3iO,P)ds^0 as P^oo, and 0<7*<7(s,P), s G [0, m], P E (0, oo). 

(2.9) 

77zen model (|1.1|) - (|1.3I) admits at least one positive equilibrium solution. 
Proof. Condition (12.81) implies 

« ( o)= fexp/- r^44 r ft w ; 8i(y ' o) expired,-) d, ds 



o 70,0) J 7 7(^,0) IJo 7^,0) 







7(2/. 0) , 

7(2,0) 

> 1. (2.10) 



Condition (12.91) and the growth behaviour of the functions in (12.71) imply that 

lim R(P) = 0, 

hence the claim holds true on the grounds of the Intermediate Value Theorem. □ 
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2.2. The general case 

For a fixed P £ (0, oo) we define the operator Bp by 

B P u = -j- ( T (., P) W ) - P)« + J /?(-, P) U (j/) dy, 

Dom(i3p) ={ue W hl (0, m) \ u{0) = 0} . (2.1 1) 

Our goal is to show that there exists a P* such that the operator Bp t has eigenvalue with a 
corresponding unique positive eigenvector. To this end, first we establish that Bp is the generator of 
a positive semigroup. Then we determine conditions that guarantee that it generates an irreducible 
semigroup. We also establish that the governing linear semigroup is eventually compact, which 
implies that the Spectral Mapping Theorem holds true for the semigroup and its generator, and the 
spectrum of the generator may contain only isolated eigenvalues of finite algebraic multiplicity (see 
e.g. [9|). It then follows that the spectral bound is a dominant (real) eigenvalue \p of geometric 
multiplicity one with a corresponding positive eigenvector flU Chapter 9]. Finally we need to 
establish conditions which imply that there exist a P + £ (0, oo) such that the spectral bound 
s(B P +) is negative and therefore the dominant eigenvalue \ P + = s(B P +) is also negative; and a 
P - £ (0, oo) such that this dominant eigenvalue X P - = s(Bp-) is positive. Then it follows from 
standard perturbation results on eigenvalues (see e.g. lfT5lO that there exists a zero eigenvalue. A 
similar strategy was employed in [2j to establish the existence and uniqueness of an equilibrium 
solution of a cyclin structured cell population model. 

Lemma 2. For every P £ (0, oo) the semigroup T{t) generated by the operator Bp is positive. 
Proof. We rewrite (12.111) as, B P = Ap + Cp, where 

d 

A P u = - — ( 7 (-, P)u) - n(-,P)u 
Dom(^lp) = {u £ W hl (0, m) \ u(0) = 0} , 

C P u= (3(-,y,P)u(y)dy, 
Jo 

Dom(Cp) =L 1 (0,m). (2.12) 
For < / £ L 1 (0, m) the solution of the resolvent equation 

(Xl-A P )u = f, 

M(s) =/ exp rX — ^) — d i^) dy - 

This shows that the resolvent operator 1Z(\, Ap) is a positive bounded operator, hence Ap gene- 
rates a positive semigroup. Since C P is a positive and bounded operator, the statement follows. □ 
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Lemma 3. The linear semigroup T(t) generated by the operator Bp is eventually compact. 

Proof. We note that Ap generates a nilpotent semigroup, while it is easily shown that Cp is a 
compact operator if conditions (11.41) hold true. (For more details see also Theorem [T2~l) □ 



Lemma 4. Assume that for every P G (0, oo) there exists an Sq > such that for all < e < £o 

/ / P(s,y,P)dyds>0. (2.13) 

JO Jm-E 

Then the linear semigroup T(t) generated by the operator Bp is irreducible. 

Proof. We only need to show that under condition (12.131) for every p G L\(0, m) there exists a t 
such that 

suppT(t )p = [0,7ft], 
for all t >t . Since 7 > 0, there exists such that 

suppT(t)p n supp /3(s, ■ ) 7^ 

for every < t and every s G (0, e]. By assumption (12.131) . T(t)p (s) > for < t and 
s G (0, e]. After this, eventually the support of the solution T(t )p will cover the entire size space 
[0,m]. □ 



Lemma 5. Assume that there exists a (3 (s, y, P) = (3 X (s, P)f3 2 {y) and a P G (0, 00) such that 

fc(a,P-)fe(y) < P(s,y,P-), s,ye [0,m], (2.14) 

and 

ffto,J-)*(.) f_ p,(,,p-)+M*,Pj d 4 iyis>l} (2 . 15) 



JO 



and a /3 + (s, P) = @i(s, and a P + G (0, 00) smc/i that 

P(s,y,P + )<Pt(s,P + )Pt(y), (2.16) 

and 



</o 



7(1/, ^ + ) I y„ 7(^^ + 



77zen ?/ze operator Bp- has a dominant real eigenvalue Ap- > and the operator Bp+ has a 
dominant real eigenvalue Xp+ < 0, with corresponding positive eigenvectors. 
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Proof. First assume that there exists a @~(s, y, P) = P)/?j(y) and a P~ such that conditions 
(12.141 ) and (12.151 ) hold true. Let £> p _ denote the operator that corresponds to the fertility j3~ and 
the constant P~. The solution of the eigenvalue problem 

Bp„u = \u, u(0) = (2.18) 



is 

(2.19) 

We multiply equation (12.191) by (3% and integrate from to m to arrive at the characteristic equation 

ft {s) l ~^r exp rl — ^ — dz \ dyds - <2 ' 20) 

Equation (12.201) admits a unique dominant real solution A p _ > if condition (12.151) holds true. 
Since £> p _ is a generator of a positive semigroup and (B P - — £> p _) is a positive (and bounded) 
operator by condition (12.141) . it follows that B P - has a dominant real eigenvalue \ P - > A p _ > 0, 
see e.g. [9, Corollary VI. 1.1 1]. 

In a similar way, let us assume that there exists a /3 + (s, y, P) = Pi(s, P)(3l$(y) and a P + such 
that condition (12.161) and (12.171) hold true. Let £> p+ denote the operator which corresponds to the 
fertility [3 + and the constant P + . The solution of the eigenvalue problem 

B% + u = Xu, u(0) = (2.21) 

is now 

u(a)= L *W"W ds x ^p^^rx — — iz \ iy - 

(2.22) 

We multiply equation (12.221) by (3% and integrate from to m to arrive at the characteristic equation 

ft (S) X K^ ew \-J, Az ) dyds - <2 ' 23) 

Equation (12.231) admits a unique dominant real solution A p+ < if condition (12.171) holds true. 
Since B P + is a generator of a positive semigroup and (£> p+ — B P +) is a positive operator by 
condition (12.161 ), it follows that B P + has a dominant real eigenvalue \ P + < A p+ < 0. 

In both cases, the positivity of the corresponding eigenvector follows from the irreducibility of 
the semigroup T[t), see 0, Theorem 9.1 1]. □ 

Theorem 6. Assume that conditions (12.131) . (I2.14I) - (I2.17I) are satisfied. Then system (|1.1I) - (11.3I) 

admits at least one positive equilibrium solution. 

Proof. Let P* > be such that s(B* P ) = 0. Then, since the spectrum consists only of isolated 
eigenvalues we have X P * = s[B P *) = and there exists a corresponding positive eigenvector p*. 
Then n—jr- p* is the desired equilibrium solution with total population size P*. □ 
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3. The linearised semigroup and its regularity 

Here, when we use the term 'linearised semigroup', we refer to the linear semigroup governing 
the linearised system. However, since it was proved in j3]| that model (ll.lh - dl.3l) is well-posed, 
there exists a semigroup of nonlinear operators £(£) t > defined via E(t)p(s, 0) = p(s, t). It was 
proven in [8 ] that if the nonlinearities are smooth enough (namely, the vital rates are differentiable) 
then this nonlinear semigroup E(i) is Frechet differentiable and the Frechet derivative around an 
equilibrium solution defines a semigroup of bounded linear operators. In this section we will 
establish the existence of this semigroup and at the same time arrive at a condition which guarantees 
that it is positive. 

Given a positive stationary solution p* of system (11.1 ft - <| 1 .3ft , we introduce the perturbation 
u = u(s, t) of p by making the ansatz p = u + p*. A Taylor series expansion of the vital rates gives 
the linearised problem (see e.g. ifTOl ) 

u t (s, t) = -7(s, P*) u s (s, t) - (j 8 (s, P*) + fi(s, P*)) u(s, t) 

- (7 s p(s,P*)p*(s) + //p(s,P*)p*(s) +7p(s,P*)p*'(s)) U(t) 

u(y,t)\B(s,y,P m ) + J (3 P (s,z,P*)p*(z)dz) dy, (3.1) 

7(0,^(0,^0 (3.2) 

where we have set 

U(t) = / u(s,t)ds. (3.3) 



Eqs. (13 . 1 b — (13 .21) are accompanied by the initial condition 

u(a,0) = «o(s). (3.4) 

Our first objective is to establish conditions which guarantee that the linearised system is governed 
by a positive semigroup. To this end, we cast the linearised system (13 . 1 1) - (I3 .4-1) in the form of an 
abstract Cauchy problem on the state space X = L 1 (0, m) as follows 

d 

— u=(A + B + C + V)u, u(0) = uq, (3.5) 

LLC 



where 



Au = -7(-, P*) u s with domain Dom(^l) = {u e W hl (0, m) | u(0) = 0} , (3.6) 

Bu = -(y,(;P.)+n(;P,))u on X, (3.7) 

pm 

Cu = - (7 s p(-, P*)p* + /i P (-, P*) p* + 7p(- 3 P*) p'J / n(s)ds 

= -p,(-) / n(s)ds on AT, (3.8) 
Jo 

pm / pm \ 

Pn= / u{y)(p{;y,P*)+ p P (-, z, P*)p*{z) dz ) dy onX, (3.9) 
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where is defined via equation (I3.8I ). Our aim is to establish that the linear operator A + B + C + V 
is a generator of a quasicontraction semigroup. To this end first we recall (see e.g. [HI SHU) some 
basic concepts from the theory of linear operators acting on Banach spaces. Let O be a linear 
operator defined on the real Banach space 3^ with norm ||.||. O is called dissipative if for every 
A > and x E Dom(O), 

\\(X- XO)x\\ > \\x\\. 
Furthermore, a function / : y — > R is called sublinear if 

f{x + y)<f(x) + f{y), x,yey 
f{Xx) = Xf(x), A>0, xey. 

If also f(x) + f(—x) > holds true for x ^ then / is called a half-norm on y. The linear 
operator O is called /-dissipative if 

f(x) < f(x- XOx), A>0, iGDom(0). 

An operator O which is p-dissipative with respect to the half norm 

p(x) = \\x + \\, 

is called dispersive, where x + = x V (and x~ = (—x) + ). Finally a C semigroup {T(t)} t>0 is 
called quasicontractive if 

||T(t)||<e^, t>0, 

for some u> E R, and it is called contractive if u < 0. We recall the following characterization 
theorem from @J. 

Theorem 7. Let y be a Banach lattice and let O : DomiO) — > y be a linear operator. Then, the 
following statements are equivalent. 

( i) O is the generator of a positive contraction semigroup. 

(ii) O is densely defined, Rg(XI — O) = y for some A > 0, and O is dispersive. 

We also recall that O is dispersive if for every x E Dom(O) there exists E y* with < </>, 
| \(ft\ | < 1 and (x, (p) = \\x + \ \ such that (Ox, 0) < 0, where (■ , ■) is the natural pairing between 
elements of y and its dual 3^*- 

Theorem 8. The operator A + B + C + V generates a positive strongly continuous ( C for short) 
quasicontraction semigroup {T(t)} t >o of bounded linear operators on X if the following condition 
holds true 

p*(s) < P(s,y,P*) +/ p P (s,y,P*)p*(y)dy, s,yE[0,m], (3.10) 
where p* is defined via equation (13.81) . 
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Proof. Our aim is to apply the previous characterization theorem for the perturbed operator 
A+B+C+V— ujX, for some w G R. To this end, for every u E Dom( v 4 + B + C + V — ujX) 
we define <ft u E X* by 



u +(s) 

~ ~ t~tt i se[0,m], u(s)^0, 



\u(s)\ 



(3.11) 



if u(s) = then let <j> u (s) = 0. Then 



and clearly 



Halloo < 1, 



u(s)4> u (s) ds = I |w H 



1- 



Making use of condition (13.101) we obtain the following estimate. 

({A + B + C + V-ul)u,<p u ) 



rm rm rm 

- l u +(s) (j(s,P*)u(s)) s ds - l u +(s)/i(s,P*)w(s)ds- / l u +(s) u; w(s) ds 
Jo Jo Jo 

rm rm / rm \ 

+ / / u(y) \P(s,y,P*)+ f3p(s,z,P s ,)p i ,(z)dz -p*(s) J dyds 



< 



l«+( s ) (t( s ) P*) u ( s )) ds - o;||m + ||i - inf /i(s, PJ lln 4 

s se[0,m] 

sup ( (3(s,y,P*) + / /3 P (s,2,P*)p*(V) dz - p*(s' 

G[0,ml V Jo 



< - (7(771, P*)«(m))l u +(m) 



sup ( /5(s, y, P*) + / 
ye[0,m] \ Jo 



/3 P (s, 2, P*)p*{z) dz - p*(s] 



< 0, 



(3.12) 



for some weR large enough, hence the operator A + B + C + V — ujX is dispersive. The operator 
A + B + C + V — 10T is clearly densely defined. We observe that the equation 

(XI-A)u = h (3.13) 

for h E X and A > sufficiently large has a unique solution u E Dom(^4), given by 



u{s) = exp 



A 



dy > / exp 



l(y,P*) J Jo {Jo l{z,P*) J i{y, p *) 



A 



dz 



Hy) 



dy. 



(3.14) 



The fact that u E Dom(^l) is well defined by (13.141) follows from 



\u'(s)\ < 



< 



h(s) 



l(s,P*) 
h(s) 



A 



+ M X , 



exp 



A 



dz 



l(z,P*) J i{y,P*) 



l%)l 



dy 
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for A large enough for some M x < oo, that is u G W 1 ' 1 ^, m). Since B + C + V — loX is bounded, 
the range condition is satisfied. Theorem[7]gives that A + B + C + V — ujX is a generator of a pos- 
itive contraction semigroup. Since the operator ujX is positive (clearly if the dispersivity estimate 
holds true with an u < then it holds true with any other u* > uj) a well-known perturbation result 
(see e.g. [9]) yields that A + B + C + V is a generator of a positive quasicontraction semigroup T 
which obeys 

\\T(t)\\<e"\ t>0. 

□ 



Remark 9. The proof of Theorem\7\shows that if 



inf fj,(s, P*) > 
se[0,m] 



sup ( /?(s,y,P*) + / (3p(s, z, P*)p*(z) dz — p*(s) j 



/zo/<i5, the growth bound u of the semigroup is negative, hence the semigroup {T(t)} t > 
is uniformly exponentially stable (see e.g. /(21/J, i.e. the equilibrium p* is locally asymptotically 
stable. 



Remark 10. We note that the operator A + B + C + V is in general a generator of a C quasicon- 
traction (but not positive) semigroup. The proof of this would utilise the Lumer-Phillips Theorem 
(see e.g. £7] @ and goes along similar lines, obtaining a dissipativity estimate in terms of u 
rather than u + , see e.g. / fiil/ . This implies that the linearised problem (I3.1I) - (I3.2I) is well-posed. 



Remark 11. Note that if (3 = /3(s, y), p — p(s), 7 = 7(5), i.e. model (ll.ll) - (ll.3l) is a linear one, 
then the biologically relevant conditions [A, (3 > and 7 > imply that it is governed by a positive 
quasicontraction semigroup. 



Theorem 12. The semigroup {T(t)} t >o generated by the operator A+B + C + V is eventually 
compact. 

Proof. C is a rank-one operator. Hence it is compact on X = L 1 (0, m). V is linear and bounded. 
Hence in view of the Frechet-Kolmogorov compactness criterion in LP we need to show that 

lim / \Vu(t + s) — Vu(s)\ ds = 0, uniformly in u, 
Jo 

for u £ B, where B is the unit sphere of L 1 (0, m). But this follows from the regularity assumptions 
we made on (3 based on the following estimate 

\Vu(si) — Vu(s2)\ < |M|i 

/3(si,y,P*) + / /3p(si, z, P*)p*(z) dz - @(s 2 ,y,P*) - / /3p(s 2 , z, P*)p*(z) dz 
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Therefore, it suffices to investigate the operator A + B. To this end, we note that the abstract 
differential equation 

^u=(A + B)u (3.15) 
at 

corresponds to the partial differential equation 

u t (s, t) + 7 (s, P m ) u s (s, t) + (7,(s, P„) + n(s, P*)) u(s, t) = 0, (3.16) 

subject to the boundary condition (13.21) . We solve easily equation (13.161) using the method of 
characteristics. For t > T(m) we arrive at 

u(s, t) = u(0, t - T{s)) exp j- f lsiy,P ^ y + ^ ,P * ] dy} = 0, (3.17) 

where 

1 



This means that the semigroup T(t) generated by A + B is nilpotent. In particular it is compact 
for t > T(m) and the claim follows. □ 

Remark 13. Theorem U2\ implies that the Spectral Mapping Theorem holds true for the semigroup 
{T{t) } t>0 with generator A + B + C + V and that the spectrum a(A + B + C + V) contains only 
isolated eigenvalues of finite multiplicity (see e.g. /(9]/J. 



4. (In) Stability 

Here, we consider the stability of positive equilibrium solutions by studying the point spectrum of 
the linearised operator A + B + C + V. The main difficulty is that the eigenvalue equation 

(A + B + C + V- J) A = 0, 

cannot be solved explicitly, since in general, the operator V has infinite rank. We encountered 
this problem previously with hierarchical size- structured population models [fTTl[T2l| . In |fTT]| and 
lfT2l we used the dissipativity approach, presented in the previous section, to establish conditions 
which guarantee that the spectral bound of the linearised semigroup is negative. However, as we 
can see from Remark [9] this approach gives a rather restrictive stability condition. Therefore, here 
we devise a different approach, which uses positive perturbation arguments. 

Theorem 14. Assume that there exists an e > such that 

P{s,y,P*) - P*(s) -e+ / Pp(s,y,P*)p*(y) dy > 0, s,ye[0,m}, (4.1) 
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and 



e / exp 
'o 



, ,-y 7s(<t,P*)+ai(c 

>.) d(T l .Jo ;W,P.) 



| j r . fa P.H : C.. f ., d<7 | 



'o 7(0-, -f*) J ^0 7(y,P*) 

TTzen the stationary solution p*(s) of model dl.lh - dl.3h is linearly unstable. 
Proof. Let e > 0, and define the operator JF £ on X as 

Jvtt = e u(s) ds = eu. 



dyds>l. (4.2) 



We first find the solution of the eigenvalue equation 

(A + B + .F b )m = Xu 



as 



= e u exp 



A + %(a, P*) + /i(a, P*) 



da 



x / — exp 

/o 712/, P) Uo 



7(^,P*) 
y A + 7 s (a,P*) +[i{a,P*) 



l&P 

Next we integrate the solution (14.3b over [0, m] to obtain 

A + 7 s (o-, P*) + /j(a, P*) 



da > dy. 



(4.3) 



u=eu 



x 



exp 



7(3/> F *) ^ t 



da 

7(<7,P*) 
y A + 7 8 (<7,P t ) + //(<r,P») 



da > dy 



ds. 



(4.4) 



We note that, if u — then equation (14.31) shows that u(s) = 0, hence we have a non-trivial 
eigenvector if and only if u 7^ and A satisfies the following characteristic equation 



A + 7s(cr, P.) + p.) 
7(<7,P*) 



da 



It is easily shown that 



1 f /• y A + 7 s (a,P,) + /J (q,P) 
7(^P) eXP U 7(^) Md?/ 



lim K(A) = 0, 

A— >+oo 



ds. (4.5) 



therefore it follows from condition (14.21) . on the grounds of the Intermediate Value Theorem, that 
equation (14.51) has a positive (real) solution. Hence we have 



< s(A + B + T e 
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Next, for a fixed < / G X, we obtain the solution of the resolvent equation 

(XL-(A + B + T E ))u = f, 

as 

u{s)=expt- — — da 

I Jo 7(<r> p *) 

x L exp u — ^ — da i dy - (4 - 6) 

We integrate equation (14.61) from to m to obtain 

Jo" exp { - j; ^gjgg g d ff } j; exp {ff ^<%}ff--*> <y}^h<iy 

It follows from the growth behaviour of the exponential function and from assumptions ( 11.41 ), that 
u is well-defined and non-negative for any < / G X and A large enough. Hence the resolvent 
operator 

TZ(X, A + B + F £ ) = (X-(A + B + F e ))~ l 
is positive, for A large enough, which implies that A + B + T e generates a positive semigroup (see 
e.g. HU). 

Finally, we note that condition (14.11) guarantees that the operator C + V — T e is positive, hence 
we have for the spectral bound (see e.g. Corollary VI. 1.11 in |[9]|) 

< s(A + B + T £ ) < s(A + B + T e + C + V-J 7 £ ) = s(A + B + C + V), 

and the result follows. □ 
Next we show that for a separable fertility function we can indeed explicitly characterize the 
point spectrum of the linearised operator. 

Theorem 15. Assume that (3 (s,y, P) = 0i(s, P)02(y), s,j/6 [0, m], P G (0, oo). Then for any 
A G C, we have A G a (A + B + C + V) if and only ifX satisfies the equation 

K,(X)=det( 1 +^ ? 2 | A) m )=0, (4.8) 
y a 3 (X) 1 + a 4 (A) J 

\ ds, 



where 



ai(A) = - 


/ F(X,s,P*) 
Jo » 


[1 


02(A) = - 


pm 

/ F(X,s,P*) 
Jo j 


L 1 


03(A) = - 


pm 

/ ft(s)F(A,s : 
Jo 


p*) 


04(A) = — 


pm 

/ &(s)F(A,s : 


p*) 



F(A,y,P 
7 (y,P,)^(A,l/ } -P.) 



dyds, (4.9) 
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and 

fo{y)p*{y) dy - p*( s ) 

» w = ^IWK) ' se[0 ' m1 ' 

F(K s, P.) = ^{-[ X + P ' ] dv] ■ .6[0,m]. 

Proof. To characterize the point spectrum of^4 + i3 + C + Pwe consider the eigenvalue problem 

{A + B + C + V - XX)U = 0, C/(0) = 0. (4.10) 
The solution of (14.101 ) is found to be 

U( 8 )=UF{\, 8 ,P.) f ° {v) p . dy + UF(\,s,P*) f Pl [ y F P * ] — dy, (4.11) 

Jo F{\y, p *) Jo i{y,P*)F{A,y,P*) 

where 

U= U(s)ds, U= (3 2 (s)U(s) ds. 
Jo Jo 

We integrate equation (14.111) from zero to m and mulitply equation (14.111) by (3 2 (s) and then inte- 
grate from zero to m to obtain 

U(l + ai (X)) + Ua 2 (X) = 0, (4.12) 
Z7a 3 (A) + U(l + 04(A)) = 0. (4.13) 

If A G er(y4 + B + C + V) then the eigenvalue equation (14.101) admits a non-trivial solution U 
hence there exists a non-zero vector (U, U) which solves equations (I4.12I) - (I4.13I) . However, if 
(U, U) is a non-zero solution of equations (I4.12I) - (I4.13I) for some A G C then (14.1 II) yields a 
non-trivial solution U. This is because the only scenario for U to vanish would yield 

^ (A > s) [ Mr) dy = - 5F(A - s) I s £ |0 - m] - 

This however, together with equations (I4.12I) - (I4.13I) would imply U = U — 0, a contradiction, 
hence the proof is completed. □ 

Theorem 16. Assume that condition (13.101) /ioWs true for some stationary solution p*. Moreover, 
assume that there exists a function j3(s, y, P) = [3i(s, P)f3 2 {y) sucn that (3(s, y, P*) < f3(s, y, P*) 
for s, y G [0, m] and characteristic equation K% (A) = Joe* not have a solution with non- 
negative real part. Then the equilibrium solution p* is linearly asymptotically stable. 
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Proof. We need to establish that the spectral bound of the linearised operator A + B + C + V is 
negative. To this end, we rewrite the operator V as a sum of two operators, namely V = Q + Hp, 
where 

pm pm 

Gu= u(y)dy (3 P (-, z, P*)p*(z) dz, on X, 
Jo Jo 

pm 

Hpu = I u(y)[3(-,y,P*)dy, on X. 
Jo 



Condition (13.101) guarantees that A + B + C + Q+Hp is a generator of a positive semigroup, while 
the eventual compactness of the linearised semigroup assures that the spectrum of A + B + C + Q+ 
Tig contains only eigenvalues and that the Spectral Mapping Theorem holds true. Since Hg — Hp 
is a positive and bounded operator we have 

s{A + B + C + g + H B ) < s{A + B + C + g + Hp + H~ p -Hp) = s(A + B + C + G + Hp) < 0, 

(4.14) 

and the proof is completed. □ 

Example 17. As we can see from equations (|4.8I) -( |4.9I) the characteristic function Kz(X) is rather 
complicated, in general. Therefore, here we only present a special case when it is straightforward 
to establish that the point spectrum of the linear operator A + B + C + Q + does not contain 
any element with non-negative real part. In particular, we make the following specific assumption 

&(•) = 

In this case we can cast the characteristic equation (14.81) in the simple form 

ex f r X + 7 ,(r, P.) + P.) d l ( g(y)i(y, P.) + ^W ] dyds = h 



JO 



7( r »P«) J V i(y,P*) 

(4.15) 

»(2/)7(2/, P*) + P*)/3 2 > 0, ye [0, m], 

which is equivalent to the positivity condition (13.101) . then equation (14.151) admits a dominant 
unique (real) solution. On the other hand, it is easily shown that this dominant eigenvalue is 
negative if 



m ps 



JO 



7.(r, P.) ± ,(r, P A | /,(,H,- P ; ) + A(y, ^ < L 

7(r, P*) J V 7 y, ^ 



(4.16) 

// is easy to see, making use of equation (12.71) . that (14.161) is satisfied if 
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holds true. In this case, we obtain for the growth bound of the semigroup u 

u = s(A + B + C + Q + H~) < 0, 

see e.g. Theorem 1.15 in Chapter VI of f^j, which implies that the equilibrium solution is linearly 
stable. 

5. Concluding remarks 

In this paper, we analysed the asymptotic behaviour of a size-structured scramble competition 
model using linear semigroup methods. We are motivated by the modelling of structured macro- 
parasites in aquaculture, specifically the population dynamics of sea lice on Atlantic salmon pop- 
ulations. First we studied existence of equilibrium solutions of our model. In the case when 
the fertility function is separable, we easily established monotonicity conditions on the vital rates 
which guarantee the existence of a steady state (Proposition [T]). In the general case we used posi- 
tive perturbation arguments to establish criteria that guarantee the existence of at least one positive 
equilibrium solution. Next, we established conditions for the existence of a positive quasicontrac- 
tion semigroup which governs the linearised problem. Then we established a further regularity 
property of the governing linear semigroup which in principle allows to study stability of equilib- 
ria via the point spectrum of its generator. In the special case of separable fertility function we 
explicitly deduced a characteristic function in equation (14.81) whose roots are the eigenvalues of 
the linearised operator. Then we formulated stability/instability results, where we used once more 
finite rank lower/upper bound estimates of the very general recruitment term. It would be also 
straightforward to formulate conditions which guarantee that the governing linear semigroup ex- 
hibits asynchronous exponential growth. However, this is not very interesting from the application 
point of view, since the linearised system is not necessarily a population equation anymore. 

Characterization of positivity using dispersivity resulted in much more relaxed conditions than 
those obtained in |[T0l for a more simple size-structured model with a single state at birth by 
characterizing positivity via the resolvent of the semigroup generator. This is probably due to the 
different recruitment terms in the two model equations. Positivity is often crucial for our stability 
studies, as was demonstrated in Section 3. Indeed, more relaxed positivity conditions result in the 
much wider applicability (i.e. for a larger set of vital rates) of our analytical stability results. 

Due to the fact that the positive cone of L 1 has an empty interior, characterizations of positivity 
such as the positive minimum principle (see e.g. HI) do not apply. However, there is an alternative 
method, namely the generalized Kato inequality (see e.g. 0]). In our setting the abstract Kato- 
inequality reads 

S u (A + B + C + V)u<(A + B + C + V)\u\, (5.1) 
for u E Dom(^l + B + C + V), where S u is the signum operator, that is 
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Inequality (15.11 ) requires 

S u J u(y) (p{s,y,P*) + J (3(s,z,P*)p*(z)dz - p*(s)j dy 

/m / rm \ 

\u(y)\\0(s,y,Pj + J 0(s,z,P*)p*(z)dz- p*(s)J dy, se[0,m], (5.2) 

which holds true for every u E Dom(*4. + B + C + V) indeed when condition (13.101) is satisfied. 

As we have seen previously in Section [XI since the linearised system is not a population model 
anymore, the governing semigroup is not positive unless some additional condition is satisfied. 
However, it was proven in |fT6l that every quasicontraction semigroup on an L 1 space has a minimal 
dominating positive semigroup, called the modulus semigroup, which itself is quasicontractive. 
Hence, in principle, one can prove stability results even in the case of a non-positive governing 
semigroup, by perturbing the semigroup generator with a positive operator such that the perturbed 
generator does indeed generate a positive semigroup. 

Acknowledgements 

JZF is thankful to the Centre de Recerca Mathematica and to the Department of Mathematics, Uni- 
versitat Autonoma de Barcelona for their hospitality while being a participant in the research pro- 
gramme "Mathematical Biology: Modelling and Differential Equations" during 01/2009-06/2009. 
PH thanks the University of Stirling for its hospitality. We also thank the Edinburgh Mathematical 
Society for financial support. 

References 

[1] W. Arendt, A. Grabosch, G. Greiner, U. Groh, H. P. Lotz, U. Moustakas, R. Nagel, 
F. Neubrander and U. Schlotterbeck, One-Parameter Semigroups of Positive Operators, 
Springer- Verlag, Berling, (1986). 

[2] R. Borges, A. Calsina and S. Cuadrado, Equilibria of a cyclin structured cell population 
model, Discrete Contin. Dyn. Syst., Ser. B 11 (2009), 613-627 '. 

[3] A. Calsina and J. Saldana, Basic theory for a class of models of hierarchically structured 
population dynamics with distributed states in the recruitment, Math. Models Methods Appl. 
Sci. 16 (2006), 1695-1722. 

[4] Ph. Clement, H. J. A. M Heijmans, S. Angenent, C.J. van Duijn, and B. de Pagter, One- 
Parameter Semigroups, North-Holland, Amsterdam 1987. 

[5] M. J. Costello, Ecology of sea lice parasitic on farmed and wild fish, Trends in Parasitol. 22 
(2006), 475-483. 



19 



J. Z. Farkas et al. 



Semigroup analysis of structured parasite populations 



[6] J. M. Cushing, An Introduction to Structured Population dynamics, SIAM, Philadelphia 
(1998). 

[7] O. Diekmann and M. Gyllenberg, Abstract delay equations inspired by population dynamics, 
in "Functional Analysis and Evolution Equations" (Eds. H. Amann, W. Arendt, M. Hieber, F. 
Neubrander, S. Nicaise and J. von Below), Birkhauser, (2007), 187-200. 

[8] O. Diekmann, Ph. Getto and M. Gyllenberg, Stability and bifurcation analysis of Volterra 
functional equations in the light of suns and stars, SIAM J. Math. Anal. 39 (2007), 1023- 
1069. 

[9] K.-J. Engel and R. Nagel, One-Parameter Semigroups for Linear Evolution Equations, 
Springer, New York 2000. 

[10] J. Z. Farkas and T. Hagen, Stability and regularity results for a size- structured population 
model, J. Math. Anal. Appl. 328 (2007), 119-136. 

[11] J. Z. Farkas and T. Hagen, Asymptotic analysis of a size-structured cannibalism model with 
infinite dimensional environmental feedback, Commun. Pure Appl. Anal. 8 (2009), 1825- 
1839. 

[12] J. Z. Farkas and T. Hagen, Hierarchical size-structured populations: The linearized semigroup 
approach, to appear in Dyn. Contin. Discrete Impuls. Syst. Ser. A Math. Anal. 

[13] P. A. Heuch and T. A. Mo, A model of salmon louse production in Norway: effects of 
increasing salmon production and public management measures, Dis. Aquat. Org. 45 (2001), 
145-152. 

[14] M. Iannelli, Mathematical Theory of Age-Structured Population Dynamics, Giardini Editori, 
Pisa (1994). 

[15] T. Kato, Perturbation Theory for Linear Operators, Springer, New York, 1966. 

[16] Y. Kubokawa, Ergodic theorems for contraction semi-groups, J. Math. Soc. Japan 27 (1975), 
184-193. 

[17] J. A. J. Metz and O. Diekmann, The Dynamics of Physiologically Structured Populations, 
Springer, Berlin, 1986. 

[18] A. G. Murray and P. A. Gillibrand, Modelling salmon lice dispersal in Loch Torridon, Scot- 
land, Marine Pollution Bulletin 53 (2006), 128-135. 

[19] J. PriiB, Stability analysis for equilibria in age-specific population dynamics, Nonlin. Anal. 
TMA 7 (1983), 1291-1313. 



20 



J. Z. Farkas et al. 



Semigroup analysis of structured parasite populations 



[20] C. W. Revie, C. Robbins, G. Gettinby, L. Kelly and J. W. Treasurer, A mathematical model 
of the growth of sea lice, Lepeophtheirus salmonis, populations on farmed Atlantic salmon, 
Salmo salar L., in Scotland and its use in the assessment of treatment strategies, J. Fish Dis. 
28 (2005), 603-613. 

[21] C. S. Tucker, R. Norman, A. Shinn, J. Bron, C. Sommerville and R. Wootten, A single cohort 
time delay model of the life-cycle of the salmon louse Lepeophtheirus salmonis on Atlantic 
salmon Salmo salar, Fish Path. 37 (2002), 107-1 18. 

[22] O. Tully and D. T. Nolan, A review of the population biology and host-parasite interactions 
of the sea louse Lepeophtheirus salmonis (Copepoda: Caligidae), Parasitology 124 (2002), 
S165-S182. 

[23] G. F. Webb, Theory of Nonlinear Age-Dependent Population Dynamics, Marcel Dekker, 
New York, (1985). 

[24] K. Yosida, Functional Analysis, Springer, Berlin, (1995). 



21 



